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Early Metabolic Markers of the Development 
of Dysglycemia and Type 2 Diabetes and Their 
Physiological Significance 
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Metabolomic screening of fasting plasma from nondiabetic 
subjects identified a-hydroxybutyrate (a-HB) and linoleoyl- 
glycerophosphocholine (L-GPC) as joint markers of insulin resis- 
tance (IR) and glucose intolerance. To test the predictivity of a-HB 
and L-GPC for incident dysglycemia, a-HB and L-GPC measure- 
ments were obtained in two observational cohorts, comprising 
1,261 nondiabetic participants from the Relationship between In- 
sulin Sensitivity and Cardiovascular Disease (RISC) study and 
2,580 from the Botnia Prospective Study, with 3-year and 9.5-year 
follow-up data, respectively. In both cohorts, a-HB was a positive 
correlate and L-GPC a negative correlate of insulin sensitivity, 
with a-HB reciprocally related to indices of (3-cell function de- 
rived from the oral glucose tolerance test (OGTT). In follow-up, 
a-HB was a positive predictor (adjusted odds ratios 1.25 [95% CI 
1.00-1.60] and 1.26 [1.07-1.48], respectively, for each standard 
deviation of predictor), and L-GPC was a negative predictor 
(0.64 [0.48-0.85] and 0.67 [0.54-0.84]) of dysglycemia (RISC) or 
type 2 diabetes (Botnia), independent of familial diabetes, sex, 
age, BMI, and fasting glucose. Corresponding areas under the 
receiver operating characteristic curve were 0.791 (RISC) and 
0.783 (Botnia), similar in accuracy when substituting a-HB and 
L-GPC with 2-h OGTT glucose concentrations. When their activity 
was examined, a-HB inhibited and L-GPC stimulated glucose- 
induced insulin release in INS-le cells. a-HB and L-GPC are 
independent predictors of worsening glucose tolerance, physio- 
logically consistent with a joint signature of IR and p-cell 
dysfunction. Diabetes 62:1730-1737, 2013 
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There is increasing interest in identifying mark- 
ers of chronic diseases. A useful biomarker 
is a molecule that 7) is easily and specifically 
measurable in accessible body fluids, 2) im- 
proves prediction algorithms, 3) tracks an underlying 
pathophysiological mechanism, 4) changes consensually 
with the mechanism, and possibly, 5) maps onto a new 
disease pathway. This search hopes not only to provide 
a more accurate prediction of disease but also follow its 
evolution and response to intervention and expose new 
potential therapeutic targets. 

Recent studies have used different extents of metab- 
olomic profiling to uncover metabolic signatures of obesity 
(1-3), fatty liver disease (4), and type 2 diabetes (T2D) (5). 
In previous work (6), we screened several hundred 
metabolites to identify novel circulating biomarkers of 
insulin resistance (IR) — as measured by the euglycemic 
clamp technique — in selected subjects in the Relationship 
between Insulin Sensitivity and Cardiovascular Dis- 
ease (RISC) study, a cohort of well-phenotyped nondia- 
betic individuals. The top biomarker of IR emerging from 
this analysis was a previously unrecognized metabolite 
identified as a-hydroxybutyrate (a-HB), an organic acid 
positioned at an interesting crossroad of intermedi- 
ary metabolism — amino acid catabolism and glutathione 
synthesis — and upstream to the tricarboxylic acid (TCA) 
cycle. The next top-ranking biomarker was linoleoyl- 
glycerophosphocholine (LrGPC), an independent correlate 
of insulin sensitivity and a putative lipid-signaling molecule. 

Because these biomarkers of insulin sensitivity were 
independently associated with glucose intolerance (6) and 
IR is a risk factor for T2D (7), we carried out a clinically 
validated assay of these metabolites in the entire RISC 
cohort at baseline and follow-up to evaluate their pre- 
dictivity of incident dysglycemia. Furthermore, because 
these biomarkers have not been previously evaluated in 
clinical outcome studies, we tested their ability to predict 
T2D in a long-term observational cohort of at-risk subjects 
in the Botnia Prospective Study. Finally, we initiated ad- 
ditional in vivo and in vitro studies to gain understanding 
of the physiological basis for these associations. 

RESEARCH DESIGN AND METHODS 

Study cohorts. The RISC study is a prospective, observational cohort study 
(n = 1,308) for which rationale and methodology have been published (8). In 
brief, participants were recruited from people attending clinic and from lab- 
oratory personnel at 19 centers in 13 countries in Europe, according to the 
following inclusion criteria: men or women, aged between 30 and 60 years 
(stratified by sex and age according to 10-year age groups), and clinically 
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healthy. Initial exclusion criteria were treatment for obesity, hypertension, 
lipid disorders or diabetes, pregnancy, cardiovascular or chronic lung disease, 
weight change of >5 kg in last 6 months, cancer (in last 5 years), and renal 
failure. Exclusion criteria after screening were arterial blood pressure >140/ 
90 mmHg, fasting plasma glucose >7.0 mmol/L, 2-h plasma glucose (on 
a standard, 75-g oral glucose tolerance test [OGTT], performed in each sub- 
ject) >11.0 mmol/L or known diabetes, total serum cholesterol >7.8 mmol/L, 
serum triglycerides >4.6 mmol/L, and electrocardiogram abnormalities. 
Baseline examinations began in June 2002 and were completed in July 2005 
and included 1,538 subjects receiving an OGTT. Of these, 1,308 also received 
a euglycemic hyperinsulinemic clamp and constituted the baseline cohort; 
cross-sectional data on this cohort have been published (8). 

All subjects of the baseline cohort were recalled 3 years later, and 1,048 
(80%) participated in the follow-up evaluation (8-10). The baseline anthro- 
pometric and metabolic characteristics of the 260 subjects who were lost to 
follow up were very similar to those of the subjects who participated (data not 
shown). The follow-up study included all the baseline anthropometric meas- 
urements and the OGTT. Fasting plasma samples were also obtained from 47 
T2D patients (14 men, 33 women; age, 56 ± 7 years; BMI, 29.9 ± 4.2 kg/m 2 ; 
HbA lc , 7.4 ± 0.8%; fasting glucose, 9.3 ± 1.9 mmol/L) after 4 weeks of oral 
hypoglycemic agents washout. 

The Botnia study is a family-based, observational familial study started in 1990 
on the West coast of Finland aiming at identifying diabetes susceptibility genes 
(11). The prospective part included 2,770 nondiabetic family members and/or 
their spouses (1,263 men and 1,507 women; mean age, 45 years), 151 of whom 
developed T2D during a median 9.5-year follow-up period (7). All subjects were 
given information about exercise and healthy diet and exposed at 2- to 3-year 
intervals to a new OGTT. Baseline measurements of a-HB and L-GPC were 
obtained from 2,580 Botnia participants, with 9.5-year median follow-up data 
analyzed. 

Ethics committee approval was obtained by recruiting centers. 
Measures. Blood samples were taken before and 30, 60, 90, and 120 min (30, 60, 
and 120 min in Botnia) after a 75-g OGTT, with a follow-up test repeated in RISC 
and Botnia. In RISC, clamps and OGTTs were performed in all subjects within 
same week (10). 

Targeted metabolite analysis. For absolute quantitation, metabolites were 
analyzed by an analytically and clinically validated isotope-dilution ultra-high- 
performance liquid chromatographic tandem mass spectroscopy (UHPLC-MS/ 
MS) assay developed and conducted in a Clinical Laboratory Improvement Act/ 
College of American Pathologists-accredited laboratory, as a consolidated 
version of the assay described previously (4). In brief, 50 jxL EDTA plasma 
samples were spiked with internal standards and subsequently subjected to 
protein precipitation by mixing with 250 fxL methanol. After centrifugation, 
aliquots of clear supernatant were injected onto an UHPLC-MS/MS system, 
consisting of a Thermo TSQ Quantum Ultra Mass Spectrometer and a Waters 
Acquity UHPLC system equipped with a column manager module in 2.5-min 
assay. a-HB, L-GPC, and oleic acid were eluted with a gradient on Waters 
Acquity single RPC-18 column (2.1 mm X 50 mm, 1.7-mm particle size) at 
a mobile-phase flow rate of 0.4 mL/min at 40°C. Ionization was achieved by 
a heated electrospray ionization source. Quantitation was performed based on 
the area ratios of analyte and internal standard peaks using a weighted linear 
least-squares regression analysis generated from fortified calibration stand- 
ards in an artificial matrix, prepared immediately before each run, with 
coefficients of variation for a-HB, L-GPC, and oleic acid at 4.0, 6.3, and 4.6%, 
respectively, based on Botnia data (n = 2,585) and 146 replicates over 9 
months. Stable deuterium-labeled compounds (a-HB-D 3 , L-GPC-D 9 , and oleic 
acid- 13 Ci 8 ) were used as internal standards. 

For targeted quantitation of amino acids, we used the Absolut/DQ pl80 kit 
(BIOCRATES Life Sciences AG, Innsbruck, Austria) on an API 4000 LC-ESI- 
MS/MS System (AB Sciex Deutschland GmbH, Darmstadt, Germany) equipped 
with an Agilent 1200 Series HPLC and an HTC PAL auto sampler (CTC Ana- 
lytics, Zwingen, Switzerland). The instrument was controlled by Analyst 1.5.1 
software. Samples were handled and preprocessed using a Hamilton Micro 
Laboratory Star robot (Hamilton Bonaduz AG, Bonaduz, Switzerland) and 
a nitrogen evaporator (Porvar, Ultravap). For data preprocessing and analysis, 
we followed the operating procedure of a proprietary method, as described in 
U.S. patent application (US 2007/0004044). In brief, 10 |xL serum were put on 
the kit 96-well plate loaded with internal standards. After derivatization of the 
amino acids with phenylisothiocyanate, metabolite extraction with organic 
extraction solvent, filtration by centrifugation, and dilution with H 2 0, the 
samples were analyzed by HPLC and MS instruments. Calculation of metab- 
olite concentrations based on internal standards was performed with the 
Met/Q package, which is an integral part of the kit. In addition, the data were 
corrected for batch effects. The Absolut/DQ methods were proven to be in 
conformance with the U.S. Food and Drug Administration Guideline (Guid- 
ance for Industry — Bioanalytical Method Validation, May 2001) requiring proof 
of reproducibility within an given error range (5). 



In vitro studies. The clonal INS-le cell line, derived and selected from the 
parental rat insulinoma, was donated by Dr. Claes Wollheim (University of 
Geneva, Switzerland). Tissue culture reagents were obtained from Gibco 
Invitrogen (Basel, Switzerland), a-HB from Sigma-Aldrich (St. Louis, MO), and 
L-GPC from Avanti Polar Lipids (Alabaster, AL). Cells were grown in mono- 
layer culture in RPMI-1640 medium containing 11.1 mmol/L glucose. The 
culture medium was supplemented with 10% heat-inactivated FCS, 1 mmol/L 
sodium pyruvate, 10 mmol/L HEPES, 2 mmol/L glutamine, 50 fxmol/L 
p-mercaptoethanol, 100 units/mL penicillin, and 100 |xg/mL streptomycin. Cells 
were cultured at 37°C in a humidified 95% air, 5% C0 2 atmosphere. Cells were 
seeded in wells at a density of 1.5 X 10 5 cells per culture well at least 96 h 
before use in the insulin secretion experiments. The response to glucose was 
tested at passage between 50 and 95. At 48 h before the experiment, when 
—80% confluence had been reached, cells were incubated in fresh RPMI-1640 
medium containing 11 mmol/L glucose, in the presence or absence of a-HB (at 
a concentration of 6, 18, or 54 |xg/mL) or L-GPC (at a concentration of 5, 15, or 
45 |xg/mL). Before experiments, cells were maintained for 1 h in glucose-free 
culture medium. Cells were washed twice and preincubated at 37° C for 1 h in 
1 mL glucose-free Krebs-Ringer bicarbonate HEPES buffer (KRBH) containing 
135 mmol/L NaCl, 3.6 mmol/L KC1, 5 mmol/L NaHC0 3 , 0.5 mmol/L NaH 2 P0 4 , 
0.5 mmol/L MgCl 2 , 1.5 mmol/L CaCl 2 , 10 mmol/L HEPES, and 0.1% BSA (glucose- 
free and free fatty acid [FFA]-free), pH 7.4. Cells were washed once with 
glucose-free KRB-HEPES buffer and subsequently incubated for 60 min with 
KRB-HEPES buffer containing 3.3, 11.0, or 20.0 mmol/L glucose, or 20 mmol/L 
glucose plus 10 mmol/L L-arginine. At the end of incubation, supernatants were 
collected to measure insulin concentration by enzyme-linked immunosorbent 
assay (Millipore Corporation, Billerica, MA). 

Data analysis. Fat-free mass (FFM) was measured by electrical bioimpedance 
(8). Glucose tolerance was categorized as normal glucose tolerance (NGT; 
fasting glucose <6.1 mmol/L and 2-h glucose <7.8 mmol/L), impaired glucose 
tolerance (IGT; fasting glucose <7.0 mmol/L and 2-h glucose >7.8 and <11.1 
mmol/L), and T2D (fasting glucose >7.0 mmol/L or 2-h glucose >11.1 mmol/L, 
or antidiabetic treatment). Impaired fasting glycemia (IFG) was defined as 
fasting glucose <7.0 and >6.1 mmol/L and 2-h glucose <7.8. IFG and IGT were 
pooled as impaired glucose regulation (IGR). 

On the basis of observed OGTT changes at follow-up, RISC participants were 
classified as stable NGT (i.e., NGT at baseline and follow-up) and progressors 
(i.e., those progressing in sequences NGT -+ IFG, NGT -+ IGT, NGT -+ T2D, 
IFG -+ IGT, IFG -+ T2D, or IGT -+ T2D between baseline and follow-up). 

Insulin sensitivity was the M value during final 40 min of clamp normalized to 
FFM ((jumol • min -1 • kg FFM _1 ). In Botnia, insulin sensitivity was estimated 
(eM) from plasma glucose and insulin concentrations measured during OGTT 
according to Stumvoll et al. (12). 

p-Cell function. In RISC, p-cell function was expressed as (3-cell glucose 
sensitivity (p-GS) as derived from mathematical modeling of plasma glucose 
and C-peptide response to glucose ingestion (13) with use of C-peptide 
deconvolution analysis from Van Cauter et al. (14). p-GS is the dose-response 
function relating insulin secretion rates to plasma glucose concentrations 
during dynamic testing. In Botnia, C-peptide measurements were not avail- 
able. Therefore, the ratio of the insulin-to-glucose (I-to-G) area under the curve 
during OGTT was used as a surrogate index of p-cell function. 
Statistical analysis. Data are presented as mean ± SD or, for variables with 
a skewed distribution, as median (interquartile range). Group values were 
compared by Mann- Whitney (or Kruskal-Wallis) test for continuous variables or 
X 2 for nominal variables, and paired group values were compared by Wilcoxon 
test. Simple associations were tested by Spearman correlation coefficient (p). 
General linear models were used to test for associations while controlling for 
potential confounders (center, sex, age, BMI), and strength of adjusted asso- 
ciations was expressed as partial correlation coefficients (r). Logistic regression 
was used to predict progression; odds ratios and 95% CIs were expressed for 
one SD change in each explanatory variable. Performance of different logistic 
models was assessed as the C statistic (area under receiver operating charac- 
teristic [ROC] curve). For in vitro data, two-way ANOVA was used to simulta- 
neously test for effect on insulin release of increasing glucose and added a-HB 
or L-GPC concentrations. Statistical analyses were run using JMP 7.0 software 
(SAS Institute, 2007). A two-sided P < 0.05 was considered significant. 

RESULTS 

Baseline. Data of RISC participants are grouped in 
Table 1 in 1,115 NGT subjects by quartile of insulin sen- 
sitivity (M value from clamp), IGR, and T2D subjects. Fa- 
milial diabetes, age, BMI, 2-h glucose, fasting insulin, and 
fasting FFA concentrations were progressively higher 
along M quartiles in NGT, IGR, and T2D subjects, whereas 
(3-GS was progressively lower. Across these groups, fasting 
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TABLE 1 

RISC study: baseline anthropometric and metabolic parameters 



NGT M quartile 





1 


2 


3 


4 


IGR 


T2D 




n = 279 


n = 280 


n = 278 


n = 278 


n = 146 


n = 47 


Worn on fO/A 


OV 


OO 


OV 


OO 


OtLt 


70* 






AA-hQ 
L ± L ± — V 


dd + 8 




d(\ + St 


DKJ — I \ 


Family history of diabetes (%)* 


18 


23 


26 


33 


43* 


80* 


BMI(kg/m 2 )$ 


23.9 ± 3.0 


24.7 ± 3.5 


25.2 ± 3.8 


27.1 ± 4.5 


27.4 ± 4.4t 


29.9 ± 4.2t 


Fasting glucose (mmol/L) 


5.02 ± 0.46 


5.01 ± 0.48 


4.99 ± 0.51 


5.00 ± 0.58 


5.49 ± 0.67t 


9.28 ± 1.89t 


Fasting insulin (pmol/L)$ 


23 (15) 


27 (17) 


30 (22) 


41 (28) 


42 (35) t 


67 (53)t 


Fasting FFA (|xmol/L)$ 


470 (240) 


470 (220) 


520 (225) 


570 (280) 


605 (310) 




M (|jimol • min" 1 • kg FFM _1 )$ 


80 (20) 


61 (13) 


49 (11) 


34 (13) 


39 (26)t 




P-GS (pmol • min" 1 • m 2 • mmol/L" 1 ) 


119 (80) 


124 (85) 


118 (88) 


114 (74) 


69 (46)t 


32 (35)t 


a-HB (|ULg/mL)t 


3.51 (2.97) 


4.11 (1.88) 


4.36 (2.12) 


4.87 (2.33) 


5.10 (2.55)t 


7.42 (2.41)t 


L-GPC (|xg/mL)$ 


16.81 (6.81) 


15.84 (5.88) 


14.90 (6.10) 


13.03 (5.18) 


13.83 (7.73)t 


7.34 (3.31)t 



Entries are mean ± SD or median (interquartile range), unless otherwise indicated. *P < 0.05 by x 2 test. fP — 0.05 for the difference between 
IGR and NGT by Mann-Whitney test. %P < 0.0001 for NGT quartiles by Kruskal-Wallis test. 



levels of a-HB were progressively higher, whereas L-GPC 
concentrations showed an inverse gradient, with re- 
spective highest and lowest levels observed in T2D. 

After adjusting for sex, age, BMI, and study site, a-HB 
and L-GPC concentrations were directly and inversely as- 
sociated with insulin sensitivity (M), respectively (adjusted 
r 2 = 0.33 and 0.34, respectively, both P < 0.0001; Supple- 
mentary Fig. 1). In the same model, a-HB was reciprocally 
related to (3-cell function ((3-GS, partial r = -0.11, P = 
0.0002), whereas L-GPC was unrelated. By defining IR as 
an M value in the bottom quartile of NGT subjects (<39 
(xmol • min -1 • kg FFM _1 ), a-HB concentration in the top 
quartile of its own distribution (>5.48 |xg/mL) confers an 
IR risk of 2.84 (95% CI 2.04-3.95), whereas an L-GPC 
concentration in the bottom quartile of its distribution 
(<11.78 |xg/mL) confers a risk of 3.14 (2.19-4.52). In the 
122 NGT subjects falling in the highest a-HB quartile and 
the lowest L-GPC quartile, the risk of IR is 4.14 (2.60-6.70). 
Prediction of IR (defined as above) increases from a ROC 
of 0.801 when using familial diabetes, sex, age, and BMI as 
predictors to a ROC of 0.837 upon adding a-HB and L-GPC 
measurements. 



Data of Botnia participants are given in Table 2 by 
quartile of eM (Stumvoll index) for the 1,811 NGT subjects 
and separately for 642 IGR subjects. By study design, 
prevalence of familial diabetes was high in this cohort and 
similar across eM quartiles in NGT. Age, BMI, fasting and 
2-h glucose concentrations, and fasting insulin concen- 
trations were progressively higher across eM quartiles and 
IGR, whereas (3-cell function (as measured by I-to-G ratio) 
was highest in the most IR quartile and lower in IGR 
subjects, indicating incipient |3-cell failure. a-HB levels 
were progressively higher across these groups, whereas 
L-GPC concentrations showed an inverse gradient. In the 
whole cohort, eM was negatively associated with a-HB and 
positively related to L-GPC (partial r = -0.19 and 0.22, re- 
spectively, both P < 0.0001) after adjusting for familial 
diabetes, sex, age, and BMI. In the same adjusted model, 
the I-to-G ratio was reciprocally related to a-HB (partial r = 
-0.09, P < 0.0001), whereas L-GPC was unrelated. 
Follow-up. In RISC (Table 3), GT was still normal in 779 
subjects (stable NGT) and had deteriorated in 123 (pro- 
gressors). The baseline clinical phenotype of progressors 
included more familial diabetes, higher age, BMI, fasting 



TABLE 2 

Botnia study: baseline anthropometric and metabolic parameters 







NGT M 


quartile 








1 


2 


3 


4 


IGR 




n = 453 


n = 453 


n = 453 


n = 452 


n = 642 


Women (%) 


57 


52 


51 


58 


52 


Age (years)* 


40 ± 13 


42 ± 12 


46 ± 13 


49 ± 12 


50 ± 13t 


Family history of diabetes (%) 


84 


82 


83 


83 


87$ 


BMI(kg/m 2 )* 


21.6 ± 2.2 


24.1 ± 1.9 


26.0 ± 2.1 


29.4 ± 3.3 


27.2 ± 4.4t 


Fasting glucose (mmol/L)* 


5.23 ± 0.48 


5.38 ± 0.46 


5.41 ± 0.44 


5.46 ± 0.42 


6.04 ± 0.56t 


2-h glucose (mmol/L)* 


4.84 ± 0.98 


5.45 ± 0.93 


5.91 ± 0.90 


6.35 ± 0.86 


7.90 ± 1.45t 


Fasting insulin (pmol/L)* 
eM (mL • min" 1 • kg" 1 )* 


30 (15) 


35 (17) 


42 (22) 


55 (35) 


52 (39)t 


10.8 (0.8) 


9.9 (0.4) 


8.9 (0.5) 


7.4 (1.5) 


7.5 (2.9)f 


I-to-G ratio (pmol/mmol)* 


4.3 (2.8) 


4.9 (2.6) 


5.7 (2.7) 


7.3 (4.5) 


5.9 (4.3)t 




3.07 (1.41) 


3.24 (1.31) 


3.41 (1.56) 


3.71 (1.49) 


3.95 (1.72)t 


L-GPC (fig/mL)* 


17.65 (5.29) 


16.61 (5.66) 


15.97 (5.03) 


14.29 (4.60) 


14.90 (5.15)t 



Entries are mean ± SD or median (interquartile range), unless otherwise indicated. *P < 0.0001 for NGT quartiles by Kruskal-Wallis test. fP < 
0.0001 for the difference between IGR and NGT by Mann-Whitney test. %P < 0.05 by x 2 test. 
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TABLE 3 

RISC study: baseline anthropometric and metabolic parameters 
by outcome 





Stable NGT 


Progressors 




n = 779 


n = 123 


Women (%) 


55 


53 




FHD (%)t 


25 


50 


Glucose (mmol/L)* 


Fasting 


5.00 ± 0.48 


5.32 ± 0.46 


2-h glucose 


5.33 ± 1.11 


6.17 ± 1.34 


Fasting insulin (pmol/L)* 


28 (20) 


39 (28) 


M (fimol • min -1 • kg FFM -1 )* 


55 (25) 


38 (28) 


p-GS 

(pmol • min -1 • m~ 2 • mnol/L -1 )* 120 (80) 69 (60) 


a-IIB (fig/mL) 


Follow-up* 


3.83 (1.73) 


4.88 (2.06) 


L-GPC (fxg/mL) 


Baseline* 


15.41 (6.60) 


13.29 (5.23) 


Follow-up* 


16.24 (7.03) 


13.02 (6.60) 



Entries are mean ± SD or median (interquartile range), unless in- 
dicated otherwise. FHD, family history of diabetes. *P < 0.0001 by 
Mann-Whitney test. fP < 0.05 by x 2 test. 



glucose, 2-h glucose, and fasting insulin concentrations. 
Furthermore, progressors were more IR and ((3-cell)- 
glucose insensitive and had higher a-HB and lower L-GPC 
concentrations. At follow-up, the a-HB had decreased 
in stable NGT subjects (by 0.27 [interquartile range, 2.00] 
(xg/mL) and increased in progressors (by 0.41 [2.1] (xg/mL), 
the difference being significant (P = 0.0003). By contrast, 
the L-GPC had increased in stable NGT (by 0.69 [6.4] |xg/mL) 
and decreased in progressors (by 0.04 [4.8] |xg/mL), this 
difference too being significant (P < 0.05; Table 3). 

Among the 2,580 Botnia participants, 151 had developed 
T2D at the 9.5-year follow-up visit (Table 4). The baseline 
clinical and metabolic characteristics of Botnia T2D pro- 
gressors versus nonprogressors were very similar to those 
of RISC progressors; again, a-HB levels were higher and 
L-GPC levels were lower. 



TABLE 4 

Botnia study: baseline anthropometric and metabolic parameters 
by outcome 





Nonprogressors 


Progressors 




n = 2,429 


n = 151 


Women (%) 


55 


53 


FHD (%) 


84 


90 


BMI (kg/m 2 )* 


25.6 ± 3.9 


28.5 ± 4.4 


Glucose (mmol/L) 


Fasting* 


5.52 ± 0.56 


5.85 ± 0.63 


2-h* 


6.15 ± 1.49 


7.45 ± 1.88 


Fasting insulin (pmol/L)* 


41 (29) 


59 (48) 


eM (mL • min" 1 • kg" 1 )* 


9.06 (2.28) 


7.18 (3.73) 


I-to-G ratio (pmol/mmol)* 


89 (95) 


71 (63) 


a-HB (fxg/mL)* 


3.49 (1.55) 


3.83 (1.62) 




15.92 (5.48) 


14.00 (3.56) 



Entries are mean ± SD or median (interquartile range), unless in- 
dicated otherwise. FHD, family history of diabetes. *P < 0.0001 by 
Mann- Whitney test. 



The predictivity of a-HB and L-GPC for incident dys- 
glycemia (RISC) or T2D (Botnia) was evaluated in multi- 
variate models including classical predictors. We found 
generally similar odds ratios between the two populations. 
BMI and fasting glucose were positive predictors in both 
cohorts, whereas familial diabetes and age were stronger 
positive predictors in Botnia than RISC. In both cohorts, 
baseline a-HB was a positive predictor and L-GPC a nega- 
tive predictor of almost superimposable strength (Fig. 1). 

We next compared the ability of a-HB and L-GPC to 
predict dysglycemia/T2D with that of traditional clinical 
models. As detailed in Table 5, adding fasting plasma glu- 
cose to the standard clinical predictors (familial diabetes, 
sex, age, and BMI) increased the ROC area under the 
curve by 0.044 in RISC and 0.017 in Botnia. Upon adding 
also the 2-h postglucose plasma glucose concentration, 
ROC area rose by a further 0.024 in RISC and 0.022 in 
Botnia. By replacing 2-h glucose with a-HB and L-GPC, 
very similar ROCs were observed in RISC (0.790) and 
Botnia (0.783). Finally, adding a-HB and L-GPC to both 
fasting and 2-h glucose improved the ROC by 0.018 in RISC 
and by 0.008 in Botnia. Thus the two biomarkers matched 
the predictivity of an OGTT (fasting insulin making a neg- 
ligible contribution) in both cohorts. Interestingly, in- 
cluding the actual measurements of insulin sensitivity and 
(3-cell function in the 2-h glucose model yielded ROC 
values of 0.817 and 0.794 in RISC and Botnia, respectively 
(i.e., only 0.016 and 0.011 higher than the fasting bio- 
markers model). 

In vitro studies. Because levels of a-HB and L-GPC were 
predictive of dysglycemia, we also tested their activity on 
insulin secretion in INS-le cells. Insulin release increased 
as glucose concentrations rose from 3.3 to 20.0 mmol/L 
and was potentiated by adding arginine to 20.0 mmol/L 
glucose. Overall, preincubation with a-HB inhibited, and 
preincubation with L-GPC potentiated, glucose- and glucose/ 
arginine-induced insulin release in a dose-dependent man- 
ner (Supplementary Fig. 2). More specifically, the effects 
of a-HB and L-GPC both appeared to be exerted on in- 
sulin secretion at low glucose concentrations. For ex- 
ample, L-GPC dose-dependently increased insulin release 
at 3.3 mmol/L glucose (from 42 ± 6 to 73 ± 13 ng/mL at the 
highest L-GPC dose, P < 0.01), and a-HB tended to dose- 
dependently decrease insulin release at 3 mmol/L glucose. 
Amino acid profile. To further examine these biomarkers 
in the context of other metabolic pathways in vivo, we 
measured amino acids and fatty acids in 542 representative 
subjects from Botnia (Supplementary Table 1), with pro- 
gressors and nonprogressors having a similar clinical 
phenotype as the entire cohort (compared with Table 2). 
Notably, branched-chain amino acids (BCAAs; leucine, 
isoleucine, valine) and three major glucogenic amino acids 
(alanine, glutamate, arginine) were increased, whereas 
glycine was significantly decreased, in progressors versus 
nonprogressors. Increased concentrations of BCAAs and 
fatty acids, such as oleate, were positively related to a-HB, 
whereas L-GPC and insulin sensitivity were reciprocally 
related to a-HB (Fig. 2). In addition, oleate and L-GPC were 
reciprocally related to one another (with partial r of -0.23, 
P < 0.0001 after adjusting for sex, age, BMI; Supplementary 
Fig. 3). 

DISCUSSION 

With T2D prevalence continuing to increase worldwide, 
there is interest in identifying high-risk patients with more 
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Sex (F) 
Familial T2D 
Age (years) 
BMI(kg/m 2 ) 
Fasting glucose (mM) 
a-HB(|jg/mL) 
L-GPC ((jg/mL) 



0.4 



~~i 1 — i — i — |~~ 

1 



Odds ratio (95% C.I.) 



FIG. 1. Multivariate logistic regression for incident dysglycemia (RISC) or T2D (Botnia). Odds ratios (95% CI) for a-HB are 1.25 (1.00-1.60) and 
1.26 (1.07-1.48), respectively, for RISC and Botnia cohorts. The corresponding odds ratios (95% CI) for L-GPC are 0.64 (0.48-0.85) and 0.67 
(0.54-0.84). Odds ratios are calculated for 1 SD of the explanatory variables. (A high-quality color representation of this figure is available in the 
online issue.) 



sensitivity than currently feasible. Because IR is a risk 
factor for IGT and T2D (15), we tested whether newly 
identified markers of insulin sensitivity and (3-cell function 
(6) could predict deterioration of glucose tolerance. We 
found a-HB and L-GPC were related to IR in an opposite 
manner — the former positive, the latter negative — in the 
RISC and Botnia cohorts, independently of known deter- 
minants of insulin sensitivity and of one another. Fasting 
concentrations of a-HB and L-GPC measured at baseline 
predicted worsening GT independently of classical pre- 
dictors and with similar power as the 2-h plasma glucose 
level. Furthermore, changes in a-HB and L-GPC over time 
tracked with their ability to predict development of 



TABLE 5 

Prediction of progression to dysglycemia (RISC study) or 
diabetes (Botnia study) 







RISC 




Botnia 


Model 


x 2 


C statistic 


x 2 


C statistic 


FHD, sex, age, BMI 


63 


0.718 


101 


0.749 


FHD, sex, age, BMI, 
fasting G 


100 


0.762 


129 


0.766 


FHD, sex, age, BMI, 
fasting G, 2-h G 


115 


0.786 


163 


0.788 


FHD, sex, age, BMI, 
fasting G, a-HB, 
L-GPC 


115 


0.790 


158 


0.783 


FHD, sex, age, BMI, 
fasting G, a-HB, 
L-GPC, 2-h G 


126 


0.804 


176 


0.796 



FHD, family history of diabetes; fasting G, fasting plasma glucose 
concentration; 2-h G, plasma glucose concentration 2 h after an 
oral glucose load. 



dysglycemia in the RISC cohort: at follow-up, a-HB had 
risen and L-GPC had fallen in subjects who progressed to 
dysglycemia compared with stable NGT subjects. 

Strikingly, odds ratios for both metabolites for risk of 
progression to dysglycemia were virtually the same in 
RISC and Botnia (Fig. 1) despite differences in population 
characteristics, sample size, and end point. Accordingly, 
the metabolites added 0.028 and 0.017 units of ROC area to 
the standard model for predicting dysglycemia in RISC and 
Botnia, respectively. In contrast, fasting insulin, which has 
been shown to predict T2D in several studies (16-18), did 
not contribute to predictivity above the joint weight of 
a-HB and L-GPC in either study. Taken together, these 
results qualify these metabolites as disease biomarkers 
through their relation to underlying pathophysiological 
mechanisms. However, their actual value in routine clini- 
cal use requires additional studies, particularly interven- 
tion studies. 

Because in RISC and Botnia an independent, inverse 
relationship also existed between a-HB levels and indices 
of (3-cell function (19), we tested the direct impact of these 
metabolites on insulin secretion. In (3-cell in vitro studies, 
a-HB dose-dependently inhibited overall glucose- and 
arginine-mediated insulin release. Although this effect 
must be confirmed in human islets and its mechanisms 
remain to be investigated, a-HB may mark IR and also 
P-cell function under diverse circumstances. This use 
would have utility for clinicians who have available sur- 
rogate indices of IR (e.g., homeostasis model assessment- 
IR [20] or Matsuda index [21]) but must choose among 
clinical tests of (3-cell function (e.g., hyperglycemic clamp, 
acute insulin response to intravenous glucose, insulinogenic 
index) that are cumbersome and relatively incon- 
gruous with one another (22-24). In this regard, it is 
notable that L-GPC appeared to exert the opposite effect 
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FIG. 2. Relationship between L-GPC, a-HB, BCAA, and oleate levels and insulin sensitivity. Fasting plasma concentrations of the BCAAs (sum of 
leucine, isoleucine, and valine) by quartile of a-HB concentrations in 542 subjects from the Botnia study are shown. Also plotted are plasma L-GPC 
and oleate concentrations, and eM values by quartile of a-HB. Plots are mean ± SEM. Note that the horizontal scale for oleate and L-GPC has been 
shifted to avoid overlapping symbols. 



to a-HB on in vitro insulin release, thereby adding to the 
evidence that certain lipid-signaling molecules (including 
lysophospholipids) can stimulate glucose-dependent in- 
sulin release (25) through lysophospholipid receptors such 
as G-protein-coupled receptor 119 (26), which localized to 
pancreatic p-cells. Interestingly, a similar lipid-signaling 
pathway has been described for G-protein-coupled re- 
ceptor 119, with its lipid agonist activation playing a role in 
glucagon-like peptide-1 release in intestinal L-cells (27). 
Moreover, lysophospholipids have been implicated in 
other signaling functions related to metabolism, including 
glucose uptake by adipocytes and myotubes (28,29). De- 
clining L-GPC levels across IR, IGR, and T2D (Table 1 and 
Table 2) are likely a reflection of an important role for this 
metabolite in glucose metabolism. 

Previous evidence has linked raised concentrations of 
selected amino acids with obesity and T2D (reviewed in 
30). Newgard et al. (1) described a metabolic signature 
related to BCAAs in obese humans. Fiehn et al. (31) 
reported increased concentrations of leucine and valine — 
in addition to lipid substrates — in a small group of obese 
T2D women. Lower glycine concentrations have been 
reported in IR offspring of T2D patients (32), consistent 
with our observation of lower glycine associated with IR 
(6). However, compared with such reported analytes, 
a-HB and L-GPC were more sensitive markers of IR, e.g., 
able to discriminate insulin sensitive from IR individuals in 
both NGT and IGR ranges (6). 

The pattern of metabolite changes we observed in the 
Botnia subgroup (Supplementarly Table 1) can be orga- 
nized into a coherent pathway (Supplementary Fig. 3) by 
assuming IR as a primary determinant of these changes. In 
this scenario, adipose tissue IR leads to elevated FFA 
concentrations (33,34), which feed into the TCA cycle and 
are oxidized at an increased rate, thereby producing an 
excess of reducing equivalents (NADH). Raised circulating 
FFA also reconstitutes phospholipids from circulating 



lipids such as L-GPC. Such overload of the TCA cycle leads 
to accumulation of amino acids such as glutamate and 
alanine, as well as a-ketobutyrate (6), substrate precursor 
of a-HB via propionyl-CoA. Oxidative stress and IR raise 
demand for glutathione synthesis (35), of which a-ketobutyrate 
and a-HB are by-products, and depletes glutathione con- 
stituents like glycine, whose levels are decreased in asso- 
ciation with IR and T2D progression (Supplementary Table 
1). Furthermore, by reducing amino acid transport (36) 
and clearance (31,37), IR also raises BCAAs (6,30,38-40), 
which also feed into the TCA cycle, directly (leucine) or 
via propionyl-CoA. Thus, increased a-HB and decreased 
L-GPC levels serve as readouts of metabolic overload 
(elevated NADH-to-NAD + ratio) and reduced glucose me- 
tabolism in both IR and the earliest phases of dysglycemia 
(6,41) (Fig. 3). Finally, the abnormal a-HB and L-GPC levels 
and their biological activity on in vitro glucose-stimulated 
insulin release may translate into a burden on in vivo p-cell 
function. 
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FIG. 3. Reconstruction of the metabolic pathway featuring a-HB and L-GPC. Unmeasured metabolites are in italic, and statistically significant 
changes between progressors and nonprogressors are indicated in red. ct-KB, a-ketobutyrate; GSH, glutathione; oc-FFA, odd-chain FFA. See text 
for further explanation. 
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